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Abstract 

We have shown recently that the resummation of soft photon contributions leads 
to a non-exponential decay of the fermion excitations in hot QED plasmas. The 
retarded propagator of a massless fermion was found to behave as Sn(t ^$> 1/gT) ~ 
exp {—aTt [lntUpt + C]}, where u> p = gT/ 3 is the plasma frequency, a = g 2 /'In, and 
C is a constant, independent of g, which was left undefined. This term is computed 
in this paper. In gauges with unphysical degrees of freedom, it is gauge-fixing 
independent provided an infrared regulator is introduced in the gauge sector. We 
also extend our analysis to hot QCD and express the quark and gluon propagators in 
the form of three-dimensional Euclidean functional integrals which may be evaluated 
on the lattice. 


*CNRS 

iCNRS 

lLaboratoire de la Direction des Sciences de la Matiere du Commissariat a l’Energie Atomique 




Figure 1: A generic n-loop diagram (here, n — 6) which is responsible for infrared di¬ 
vergences in perturbation theory. All the photon lines are soft and dressed by the hard 
thermal loop. The fermion line is hard and nearly on-shell. 

1 Introduction 


The Bloch-Nordsieck (BN) approximation [11] offers an economical description of the non- 
perturbative interactions between charged particles and soft photons. At zero tempera¬ 
ture, it provides the correct structure of the fermion propagator near the mass-shell 0. 
At finite temperature, the Bloch-Nordsieck approximation has been used, by Weldon, to 
verify the cancellation of the infrared divergences in the production rate for soft photons 
0. The remarkable structure of the “hard thermal loops” (HTL) Q emerges from similar 
kinematical approximations, as clearly emphasized in the kinetic derivation of the HTL’s 
B §■ More recently, a similar approximation has been used in Refs. [7j, Q to eliminate 
the infrared divergences in the computation of the fermion damping rate || 

As shown in Ref. ||, this calculation requires the resummation of an infinite class of 
multi-loop Feynman graphs of the type shown in Fig. |l]. These are the same diagrams as 
those of the quenched approximation (i.e., all fermion loops are ignored), except for the 
fact that the photon lines include the hard thermal loop correction. 


10 , 11 , 12 


Throughout this work, we shall be mainly interested in the leading contribution of 
such diagrams to the propagator of a hard fermion (p ^ T) near its mass-shell (p 0 ~ p)- 
As shown in Ref. §, this leading contribution can be estimated in the Bloch-Nordsieck 
approximation, that is, with the following simplified Feynman rules: (i) the fermion prop¬ 
agator 


Go(p-q) 


1 

(p° - 9°) - v • (p - q) ’ 


( 1 . 1 ) 


(ii) the photon-fermion vertex and (iii) the HTL photon propagator *D /lu (q). 

Here, = (p 0 ,p) is the external hard momentum, with p 0 ~ p, v = p/p is the corre¬ 
sponding velocity, = (1, v), and q M = (g 0 , q) is a linear combination of the soft momenta 
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of the internal photons. The above Feynman rules govern the interactions between any 
kind of hard charged (or colored) particles, irrespective of their spin, and soft gauge fields, 
to leading order in an expansion in powers of the soft momenta [T3|] (see also Appendix C 
below for a discussion of the non-Abelian case). They lead to simplifications by ignoring 
those degrees of freedom — in this case, spin and negative-energy states — which play 
no dynamical role in the kinematical regime of interest. 


The imaginary part of the fermion self-energy computed with these rules exhibits 
infrared divergences near the mass-shell, to all orders in perturbation theory. For instance, 
in the one-loop approximation we have: 


Im£o' ) (a;~p) ~ —aT In ——^—-, (1.2) 

\uj-p\ 

where a = (? 2 /47t, u> p = gT /3 (the plasma frequency), and the approximate equality 
means that only the singular term has been preserved. For two or more photon loops, the 
mass-shell divergences are power-like [§. Such divergences prevent us from computing 
the mass-shell structure of the charged particles, and in particular from obtaining the 
fermion lifetime in perturbation theory. 

Note, however, that no infrared divergences are encountered when the perturbation 
theory is carried out directly in the time representation: the inverse of the time acts then 
effectively as an infrared cutoff. For instance, the one-loop correction to the retarded 
propagator Rr(£, p) at large times is given by: 


5Sn\t, p) ~ —it [ d t'e ipt Y,ft\t',p). (1.3) 

Jo 

This expression is well defined although the limit t —» oo of the integral over t’ (which is 
precisely the on-shell self-energy £^(u; = p )) does not exist. We actually have || 


EnVt, p) ~ —iaT - for t —, 

R K ’ t UJ p 

and therefore 

5Sft\t, p) ~ —aTt I — = —aTtln(uj p t). 

J 1 j ujp C 


(1.4) 


(1.5) 


As shown in Refs. im this correction exponentiates in an all-order calculation: 


S’n(t, p) oc exp (—aTt lnuj p t 


for t 


UJrn 


( 1 . 6 ) 


Note, however, that the approximations used in Refs. Bi- and which lead to 


eq. (|1.6| ), are reliable only for computing the leading large-time behaviour displayed in 
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eq. ( |1.6| ). (These approximations involve, aside from the Bloch-Nordsieck approximation, 
also a restriction to the static photon mode.) The subleading term — i.e., the constant 
term under the logarithm — could not be obtained in this way, and the issue of its gauge 
(in)dependence remained an entirely open problem. 

In this paper, we improve the accuracy of our previous calculation by also including 
the non-static photon modes within the Bloch-Nordsieck calculation. This is sufficient 
to fix the term of order g 2 T in the exponential (F6). As we shall see, this term, which 
receives contributions from both the electric and the magnetic sectors, becomes gauge¬ 
fixing independent when an infrared regulator is introduced in the gauge sector. The final 
result, to be derived below, is 


Sn(t^> l/tv p ) oc exp j— aTt(\n(uj p t) + 0.12652... + 0(g)) j. 


(1.7) 


This result applies to a massless fermion with momentum p ~ T or larger. The extension 
to a massive (m T) test particle is straightforward. The case of a soft fermion (p gT), 
on the other hand, requires the full machinery of the HTL-resummation [[J], and will be 
not addressed here (see Ref. || for the leading order result in this case). 

In order to derive eq. (0, we shall use a finite-temperature extension of the Bloch- 
Nordsieck (BN) model, to be introduced in section 2. Formally, our construction is a 
straightforward generalization of the corresponding model at zero-temperature, as de¬ 
scribed for example in Ref. |2J. However, unlike what happens at zero-temperature, at 
finite-temperature, the BN model cannot be solved in closed form (see also [§]). The 
technical difficulty comes from the thermal boundary conditions to be imposed on the 
BN propagator, and more specifically from the thermal occupation factors for the hard 
fermion. However, as it will be explained in Sections 3 and 4, this problem can be over¬ 
come, within the desired accuracy, and this eventually yields the large-time behaviour 
indicated in eq. 0. The independence of this result with respect to the choice of the 
gauge is further analyzed in section 5. Finally, in section 6, we consider an extension of 
the thermal BN model to QCD. Because of the mutual interactions of the soft gluons, the 
non-Abelian model cannot be solved analytically. Our main result here is an expression 
of the retarded propagator of a hard quark or gluon in the form of a functional integral 
over three-dimensional Euclidean gauge fields. This representation, which is reminiscent 
of the dimensional reduction sometimes performed in the computation of static thermal 
correlation functions [b|, ITS], is well adapted to numerical calculations on a lattice. 
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Figure 2: Complex-time contour for the evaluation of the thermal expectation values: 
C = C + U C- U Co. On C. |_, t takes all the real values between t 0 to tf (eventually, we 
let tf —> cx)). On C_, t — > t — ?'0 + , where t runs backward from tf to t 0 . Finally, on C 0 , 
t = t 0 — ir, with 0 < r < (3. 

2 The Bloch-Nordsieck propagator 

As mentioned in the Introduction, we are interested in the large time decay of the prop¬ 
agator of a hard fermion moving through a QED plasma at very high temperature T : 
T 3> m e , where m e is the electron mass in the vacuum. This fermion can be either a 
thermal electron, with typical momentum p ~ T and ultrarelativistic dispersion relation 
E p = p (we neglect the electron mass relative to T), or a (generally massive) test charged 
particle, with three-momentum p ^ T and dispersion relation E p = \Jp 2 + m 2 . By test 
particle, we mean a particle which is distinguishable from the plasma particles, and is 
therefore not part of the thermal bath. The general formalism below will be developed 
for a thermal particle. We shall indicate later how one can derive from it the simpler case 
of the test particle. Also, we shall write the general formulae for a massless fermion. The 
corresponding formulae for a massive test particle will be presented only briefly. 

We are eventually interested in the retarded propagator, 

iS R {x - y) = 6(x o - y 0 ) {{ip(x ), $(y)}), (2.1) 

where the curly braces denote the anticommutator of the fermion field operators, and 
the angular braces, the thermal expectation value. However, to calculate this propagator 
at finite temperature, it is convenient to consider first the time-ordered (or Feynman) 
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propagator, 


iS(x - y) = (T i}{x)il}{y)) = 9(x 0 - y 0 ) S > (x - y) - 9{y 0 - x Q ) S < (x - y ), (2.2) 

and to observe that the two two-point functions 

S > (x-y) = {i/j(x)ijj(y )), S < (x-y) = {ip(y)i/j(x )), (2.3) 


are analytic functions of their time arguments. These functions can be first computed in 


the imaginary-time formalism |T6[, and then continued to the real-time axis. Then, the 
retarded propagator (|2.1|) can be obtained as 


iS R (x - y) = 0(x o - y 0 ) (S’fi - j /) + S < (x- j/)) . 


(2.4) 


This is a method that we have used in Ref. [[|, and we shall use it again in Appendix A 
where we compute the propagator of a hard thermal fermion. 

An alternative formalism, which permits a direct evaluation of real-time Green’s 
functions, is based on the use of an oriented contour C in the complex time plane, as 
shown in Fig. |2] [JT6[ . We define the contour-ordered propagator 


iS(x — y) = (T c ^{x)^{y)) = 9 c (x 0 , y 0 ) S > (x - y) - 9 c (y 0 , x 0 ) S < (x - y), (2.5) 


where the time variables Xo and yo lie on G, and T c and 9q denote respectively the 
contour-ordering operator and the contour theta function. (If one gives a parametric 
representation of the path, t = z(u), with u real and monotonically increasing, then path 
ordering corresponds to the ordering in u, and 9c(ti,t 2 ) = 9{u\ — u 2 ).) The contour 
propagator ( p.5|) satisfies the Kubo-Martin-Schwinger (KMS) boundary condition 

S(t 0 - y Q ) = -S(t 0 -y 0 - ifi). (2.6) 


16 


(In this equation, and often below, we omit the spatial coordinates, for simplicity.) It can 
be given the following spectral representation 0: 

d P p-ip-(x-y) 


iS(x — y) = 


( 277 ) 


'pip) 9 c (x 0 ,yo) - n(p Q ) 


(2.7) 


where p(p) is the fermion spectral density and n(p 0 ) = l/(exp(/?p 0 ) + !)• Note that, once 
the spectral density is known, the retarded propagator (|2.1|) can be obtained as: 

d P p-ip-ix-y) 


iS R (x - y) = 9(x 0 - y 0 ) 


( 277 ) 


Pip) • 


( 2 . 8 ) 


In fact, in cases where we shall use the contour method below, the relation between S 
and S R will be even simpler. Indeed, in these cases — that of a test particle, and that of 
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a very energetic thermal particle, with p T — the statistical factor n(po) can simply 
be ignored, so that eqs. (|2.7| ) and (|2.8|) become identical for x 0 and yo real. 

The large time behaviour of the fermion propagator is governed by the interactions 
of the fermion with soft thermal photons. These can be analyzed in the BN approximation. 
The propagator has then the following functional integral representation ||] (see also Ref. 

i) 


S(x-y ) = Z~ l J VAG(x,y\A)e iSc W, 
where G{x,y\A) is the solution of the equation: 

i(v • D x )G(x, y\A) = Sc(x, y), 


(2.9) 


( 2 . 10 ) 


where D M = d M + igA^, v M = (l,v), |v| = 1 , 5c(x,y) is the contour delta function 
16|], and S'cfR] is the effective action for soft photons in the hard thermal loop (HTL) 


approximation 14 hi. H 


Sc[A] = S c A ‘ x {-\ F ^F'‘“-^(d-Af^ + j c d‘xJ c d a y ^A’‘(x)n„ x (x,y)A‘'(y) 
= f c dix f c d ‘y \ AP (x)D-px-y)A , ’(y). (2.11) 


We have written this equation in the covariant gauge with parameter A. The Coulomb 
gauge V • A = 0 will also be used in what follows. is the photon polarization tensor 
in the HTL-approximation |T7J [I| . 

The gauge fields to be integrated over in eq. ( [2.9D satisfy the periodicity condition 
A M (i 0 ,x) = A^(t 0 — i/3,x.). Correspondingly, the photon contour propagator satisfies 


yo) D^ito yo ^ 4 ), 


( 2 . 12 ) 


and can be given the following spectral representation: 


-G fox ! x 


(x-y) = -i /4p e ”V^(«)[9cUo,!/o) + N(qo) 


(2.13) 


where p^ u (g) is the photon spectral density in the HTL approximation ]18| (we follow here 
the notations in Ref. [§]), and N(q 0 ) = 1/(exp(/4?o) — 1) is the Bose-Einstein statistical 
factor. 


Eq. ( |2.10|) dehnes the BN propagator of a charged particle in a classical background 
held A^x). The vector v in this equation is to be identified with the particle velocity. 
Then, the Feynman rules generated by the functional integral (|2.9| ) coincide with those 
given in the introduction. Since, in the BN model, v is a fixed parameter, the underlying 
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physical approximation is the lack of fermion recoil. This approximation is justified when 
the momentum transfer from the background field is small, at most of order gT. (The 
energy-momentum scale gT is set by the polarisation tensor 11^; see section 3 below.) 


In order to satisfy the KMS condition (|2.(j|), the BN equation ( J2.10| ) has to be solved 
for antiperiodic boundary conditions: 


G(t 0 , yo\A) — — G(t 0 — i(3 , y^A), 


(2.14) 


and similarly for y 0 . These conditions complicate the resolution of the thermal BN model 
in general. There is a simple case, however, where this complication is absent, namely 
the case of a test particle. For this case, the thermal BN model can be exactly solved, as 
we discuss now. 


2.1 The test particle 


The propagator of a test particle has only one analytic component, namely S > (x — y ); 
S K vanishes since the thermal bath acts like the vacuum for the field operators of the test 
particle. Therefore (cf. eqs. (pT|)-(p~4|)), S(x — y) — Sr(x — y) — —id(x 0 — y 0 )S > (x — y) 
and the KMS conditions (|2.6| ) do not apply. In the BN approximation, S(x,y) is still 
given by eq. but now G(x, y\A) obeys retarded conditions in real time, and not the 

conditions (|2.14| ). The solution to eq. ( |2.10| ) is then obtained in closed form 0 (see also 
Ref. HI) : 


G R (x,y\A) = -iO(x° - y°) S {3) (: 


x - y - v(x° - r/ 0 )) U(x,y) 

roc 

—i / d tS^ 4 \x — y — vt)U(x,x — vt), 

Jo 


(2.15) 


where the time variables Xo and yo are real. The three-dimensional delta function describes 
straightline propagation with velocity v. The background gauge field only contributes a 
phase factor: 


U(x, x — vt) = expj— ig j dsu ■ A(x — v(t — s))j . 

In momentum space, the free (A^ = 0) retarded propagator reads 

G£(o;,p) = - ■ , 

uj — v • p + irj 

corresponding to the following free spectral density: 


(2.16) 


(2.17) 


Po(uj, p) = ” 2ImG^(u;, p) = 2’k5{uj — v • p). 


(2.18) 
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According to these equations, the fermion mass-shell in the BN model corresponds to 
u = E p = v ■ p. Since v is to be identified with the velocity p/p of the massless fermion, 
the free mass-shell is at = p, as it should. 


In order to perform the functional integration ( [2.9P , we first rewrite the parallel 
transporter ( 2.161 ) as follows: 



L'(x.x-vt) = exp d 4 ^ jp(z)A^(z )|, 



rt 

jn{z) = gi’f, dsd (4 \z - x + v{t - s)). 

Jo 

(2.19) 

Then, 

a straightforward calculation yields: 



S R (t, P ) = -i6(t)e- u ^ A(t), 

(2.20) 

with 

A(t) = ex p| ~\f c d 4 x J c d 4 y^( x ) D ^( x -y)f(y) j, 

(2.21) 


where the time integrals run on C + only, in accordance with eq. ( 2.19 ). 

By using eqs. (|2.19|) and (|2.13|), and after a simple calculation, we can rewrite A (t) 


m 


eq. Q2.21|) as follows: 


m = exp - ^ 


d A q 


(27r) 4 Jo Jo 


d Sl [ ds 2 e- i ^ s '-^p(q)\9(s 1 -s 2 )+N(q 0 )} , ( 2 . 22 ) 


with p(q) = v^ l p fiu (q)v l/ . Note that eq. (|2.22|) could have been obtained from the cor¬ 
responding expression in the vacuum (see Ref. [J2|) by simply replacing in the latter 
the bare photon propagator by the corresponding thermal propagator for a soft photon. 
Converserly, the zero-temperature BN propagator can be obtained from eq. Q2.22| ) by 
substituting N(q 0 ) —> —0(—q o ) and replacing p^ with the free photon spectral density. 

After performing the Si and S 2 integrations, and also using the parity property 
p(—q) = —p(q), we finally cast eq. (|2.22|) into the form : 


A(t) = exp{if$(f)} A(t), 


where 


= 9 2 / 


d 4 q p(q) 


and 


(27r) 4 2(v ■ q) 
d 4 g 


A(f) = exp j-£ 2 J p{q) N(q 0 ) 


sin t(v ■ q) 
t(v ■ q) 

1 — cos t(v ■ q) 
(v ■ q) 2 


(2.23) 


(2.24) 


(2.25) 































In the derivation of the above formulae, eqs. ( |2.22| )-( P.25| ), there was no explicit 
restriction on the photon momenta g M . Since the BN model can only be trusted for soft 
photons, we need to verify that the large time behaviour of the fermion propagator, as 
given by eqs. (|2.20|)-(p.25|), is indeed controlled by soft momenta, q « T. 


In fact, the momentum integrals in eqs. (|2.24f) - (|2.25|) contain ultraviolet divergences 
coming from their zero-temperature contributions. There is a linear UV divergence in 
<&(£), and a logarithmic divergence in A (t). These divergences can be absorbed respec¬ 
tively, by mass and field-strength renormalizations |2|. However, the finite part of the 
phase <&(£) remains dominated by hard momenta contributions, and therefore is not con¬ 
sistently determined by the present approximation. Since <3>(f) does not enter the calcu¬ 
lation of the lifetime, we shall ignore it in what follows. 


The damping effects are entirely described by the function A (t), eq. (|2.25|) , which 
extends our previous result @•1 by including the effects of the non static (go 7^ 0) electric 
and magnetic field fluctuations. At T = 0, A r=0 oc exp(— g 2 In At), where A is the upper 
momentum cutoff |J. After UV renormalization, the cutoff A is replaced by the physical 
electron mass, thus yielding 


\S(t)\ ~ (mt) (3 A ' ) 2?r , 


(2.26) 


in the covariant gauge with gauge-fixing parameter A (see eq. (pi.lip ). Such a gauge- 
dependent, polynomial dependence on time merely reflects the renormalization of the 
wave-function of the fermion due to its coupling to soft virtual photons. This is not to 
be interpreted as a damping phenomenon. The mechanism which takes place at high 
temperature and which eventually gives rise to the damping of the test particle excitation 
is the exchange of soft photons between the test fermion and the thermal charged particles. 
We shall verify in section 3 that, for sufficiently large times t 1 /gT, such a collision 
involves dominantly soft photon momenta q < gT. 

For a fixed large time f, the function 


f(t,v-q ) = 


1 — cos t(v ■ q) 


(v ■ qf 


(2.27) 


in eq. (|2.25|) is strongly peaked around v ■ q = go — v ■ q = 0, with a width ~ 1/t. In the 
limit t —> oo, f(t,v ■ q) —> tt tS(v ■ g). In the absence of infrared complications, we could 
use this limit to obtain the large time behaviour of A (t). This procedure would then yield 
A (t —> oo) ~ e~ 7t , with: 


7 = ng 2 


d 4 g 

(2 


p(q) N(qo) S(v ■ q). 


(2.28) 
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We recognize in eq. ( |2.28|) the one-loop damping rate 7 = - Im (w = p ) B 


which we know, however, to be infrared divergent (cf. eq. ( |1.2|) ). Thus, in studying the 
large-time behaviour of eq. (|2.25| ), one should keep the time finite when performing the 
momentum integral. As already mentioned after eq. the inverse time plays the role 

of an infrared cutoff. This will become explicit in section 3 below. 

Eqs. ( [2.20|) -( [2.25 ) generalize trivially to a test particle with mass m. The mass-shell 
is shifted to E p = v • p + m(l — v 2 ) 1 / 2 , (which, since v = p /E p , corresponds indeed to 
E p = v^p 2 + m 2 ), and the retarded propagator has the form (|2.20| ): 


|Sfl(t,p)| = 0(t) A„(t), 


(2.29) 


where A v (t) is the function (|2.25|) with, however, |v| < 1. 


2.2 The thermal fermion 


The case of a thermal electron with momentum p ~ T is physically more interesting, since 
this is a typical quasiparticle of the plasma. Technically, however, this is more involved, 
since the KMS boundary conditions (|2.f4 ) must be taken into account. 

To appreciate the difficulty, consider the free contour propagator, as obtained by 


replacing p(p) with po{p) = 2'k5{uj — v ■ p) in eq. (| 2 . 7|) : 


G 0 (t-t',p) = 




O c (t,t')(l-n p ) - 9 c (t',t)n f 


(2.30) 


where n p = n(v • p). By using this propagator, we can solve the BN equation ( j2.1G| ) as 
a series in powers of gA p . To this aim, one can first transform eq. ( | 2 . 10 | ) into an integral 
equation, 

G(x,y\A) = G 0 (x — y) + g [ d 4 z G 0 (x - z) v ■ A(z) G(z, y\A ). (2.31) 

Jc 

Then, by iteratively solving this equation, one generates the perturbation series for 


G(x,y\A). However, in contrast to what happens for the retarded propagator (|2.15|) , 
the resulting series for the contour propagator G(x,y\A) does not exponentiate 0. The 
exponentiation of the perturbative series for Gn{x,y\A) is related to the fact that the 
retarded free propagator, 


Gg(t,p) = 


(2.32) 


satisfies the simple multiplication law G^(t, pi) G^(t, p 2 ) = — G^(t, Pi+P 2 )- The contour 
propagator Go(t, p) does not enjoy this property, because of the presence of the statistical 
factors in eq. (|2.30|) . 
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This argument suggests that the contour BN propagator may exponentiate whenever 
the fermion occupation numbers play no dynamical role. This is what happened for the 
test particle in the previous subsection, and, more generally, this will also happen for 
a thermalized fermion with very high momentum, p > T, whose thermal occupation 
number is exponentially small: n{p) ~ e -/3p 1. In fact, when n p —> 0, the free contour 

propagator ( |2.30|) reduces to the retarded function ( |2.32| ) (for real time variables). It is 
then easy to verify that the previous solution of the BN model, as given by eqs. ( [2.20| )- 
Q2.25 ), also applies to such a very energetic thermal particle, up to corrections which are 
exponentially small when p T. In particular, the case of the test particle is formally 
recovered as the limit p/T —> oo. 


What is less obvious is that the same solution holds also for a typical thermal 
fermion, with momentum p ~ T. More precisely, as will be verified in Appendix A, the 
thermal fermion propagator decays according to the same law as above, that is, 


|<S*(t,p)| oc A (t), (2.33) 

(with A (t) as defined in eq. (|2.25|) ), up to corrections of order q/T ^ g. Physically, this 
reflects the fact (which has been already mentioned at several places, and will be verified 
in the next section) that the fermion decay at large times, t 1 /gT, is determined 

by its interactions with soft photons, with momenta q ^ gT. Such interactions do not 
significantly change the electron momentum, so that the associated thermal occupation 
factors play no dynamical role. 


3 Large time behaviour 

We are now in a position to study the large-time behaviour of the fermion propagator, as 
described by the function A (t), eq. fl2.25|) . We shall verify below that the relevant energy 
scale is hidden in the photon spectral density p pv (q), and is of the order gT. Therefore, 
“large times” means times larger than 1/ gT. 

For the computation below, we shall use the photon spectral density in the Coulomb 
gauge: 

p(q o,q) = Pi{qo,q) + (l-(v-q) 2 ) pt(qo,q)- (3.1) 

(The issue of the gauge dependence will be addressed in the next section.) The two pieces 
Pi(qo,q ) and pt(qo,q) of the spectral density correspond respectively to longitudinal and 
transverse photons, which are renormalized differently by plasma effects [6]. |16j]. They 
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have the following structure (with s — l or t): 


Ps(qo , q) = 2vre(g 0 ) z a (q)5($ - u 2 s (q)) + (3 s {q 0 , q)9(q 2 - q 2 0 ), (3.2) 


and involve delta functions associated to plasma waves at time-like momenta (q$ = 
u> 2 (q) > q 2 ), and smooth contributions f3i and (3 t at ql < q 2 arising from Landau damping. 
For given go and q, the energy-momentum scale in eq. ( 13 . 2D is set by the plasma frequency 
uj p = gT/ 3, for both the on-shell and the off-shell spectral densities (see [f[8|, |], [H| for 
more details). 


For generic times, both pieces in eq. f fh2[ ) contribute to eq. 


i) That involving h(gp — u> 2 (q)) describes the emission or the absorption of an on-shell 
plasmon. By kinematics, this is only possible if the fermion is sufficiently off-shell, \u— p\ ^ 
gT: indeed, the plasmons propagate as massive particles, with (momentum dependent) 
thermal masses of order gT [17, ^ p~6| . 

ii) The contributions involving (3i and f3 t describe collisional damping, where the fermion 
exchanges a virtual photon with the other charged particles of the plasma. Such processes 
have no kinematical restrictions, and they are the only one to contribute at very large 
times t 1 /gT. 


To study the large-time behaviour, we restrict therefore ourselves to collisional pro¬ 
cesses, i.e., retain only f3i and /3 t in the photon spectral functions. From perturbation 
theory, we know that the infrared complications are related to the singular behaviour of 
the magnetic spectral density as q 0 q —► 0 || : 


— Pt{qo < q) ^ 
qo 


3n uj 2 q 

2 q 6 + (37ro;2g 0 /4) 2 


% i(®) 

qz 


as q 


0 . 


To isolate this singular behaviour, we write 


(3.3) 


— Pt(qo,q) 
qo 


2tt 6(q 0 ) 



q 2 + < 


+ —Vt(qo,q)- 
qo 


(3.4) 


A contribution oc l/(q 2 +u 2 ) has been subtracted from the singular piece — and implicitly 
included in zy(g 0 , q) — to avoid spurious ultraviolet divergences: written as they stand, 
both terms in the r.li.s. of eq. ( |3.4| ) give UV-hnite contributions. Note that by neglecting 
the regular piece zy(g 0 , q) in the right hand side of eq. O , one recovers our previous 
result in Refs. |j (as also expressed in eqs. ( |3^ ) and ( |3.9| ) below). 

With ((3.4|), the integral in eq. (|2.25|) may be separated into two pieces: 



cos 2 9 fl t (qo, q) + "t(qo,q 


cos t{v ■ q) 

(v ■ q ) 2 


(3.5) 
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and 


F m {t) = 


d 3 q ( 1 

(2?l^ 


1 — cos t(v ■ q) 

(v ■ q) 2 


(3.6) 


q z +u;J 

The first piece, F reg (t ), is infrared safe, and its large-time limit can be taken by replacing 
f(t,v ■ q ) by 7 TtS(v ■ q) (see eq. ( |2.27|) ). This yields 

Fre 9 ^ = ^J 0 dqq J_ q ^^ ^A(9o,9) - ^Pt(qo,q) + vt{qo,q)\ , (3-7) 

where we have used the delta function 5(qo — q cos 9) to perform the angular integration. 
The remaining integral occurs also in Ref. [pd]], as part of the one-loop damping rate, 
and was computed there by using sum rules plus numerical integration. It is computed 
analytically in Appendix B, with the result: 


?(t) = — In3 ~ — x 0.54931. 

87T 47T 


(3.8) 


Note that this result comes entirely from the electric piece f3i(qo,q): the two magnetic 
pieces, ( ql/q 2 ) /3t(qo, q) and u t (q 0 , q), happen to cancel each other in the final result. This 
is purely accidental, consequence of our specific choice for the substracted term l/(q 2 +u 2 ) 
in eq. fl3~3]). 

The second piece, F IR (t ), contains the potentially singular magnetic contribution, 
so that we should take the large time limit only after performing the integral over q. This 
has been done in Ref. ||, with the following result ( 7 # is the Euler constant): 

Fm{t) = j- (in<V + (7 e ~ 1) + 0(l/u p t)y (3.9) 


> 


Note that the energy scale lu p inside the logarithm arises from the large momentum (q 
gT) behaviour, where the substracted term l/(q 2 + u 2 ) acts effectively as an UV-cutoff. 

The final result for the large-time propagator reads then: 

A(t 3> l/u;p) — exp j— aTt(ln(u>pt) + 0.12652... + 0(g, l/u/pf))}. (3.10) 

For the consistency of our approximations, it is important to observe that this result has 
been obtained by integrating, in eqs. (|3.7| ) and ( fhfj ). over soft photon momenta q ^ gT. 
While this is obvious for eq. fl3.6|) , where the two terms inside the parantheses mutually 
cancel as q ^ gT , it can be also verified for eq. (B, by using the known behaviour of (3i 


and f3t at large photon momenta 18 


For completness, let us also give the corresponding results for a massive test particle: 
After reinserting the appropriate factors of v = |v| in the previous results, we get 

A v (vt 3> 1/ujp) ~ exp {—aTvt (in (uj p vt) + ( 7 # — 1) + C(u))j, (3.11) 
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where C(v) is given by the following integral: 


C(v) = 4 / / 

v 2 Jo J- 


ir^-(pi(Qo,q) - % Pt{qo, q) + v 2 v t (qo, q) 

-vq 2nq 0 \ q 2 


(3.12) 


For a very heavy particle m ^ T, we may consider the non-relativistic limit v -C 1 (this 
is consistent with our approximations as long as p ~ rnv gT). At small v, the leading 
contribution to eq. (|3.12|) comes from the electric sector. The magnetic contribution 
involves a supplementary factor of u 2 , and vanishes as v —> 0. However, because of its 
infrared sensitivity, the contribution of the magnetic sector is not analytic in v. We 
evaluate this contribution in Appendix B, where we find: 


vC(v) = i(l + i.ln^ + | + 0(ti 2 )), 


(3.13) 


where the first term, independent of v, is the contribution of the electric sector. Together 


with eq. (|3.11|) , this yields 

. , . . f aT 

A V (vt > 1 /Up) ~ exp |- —t 


1 + v In ( 4 “ v 3 (cu p t ) 2 ) + v(2'y i 


3/2) 


(3.14) 


for v 1. In particular, as v —> 0 (i.e., m —> cxj), the damping is purely exponential, 
with a damping rate 70 = aT/2 which coincides with the one-loop result in Ref. ||. As 
for the u-dependent terms, the coefficient of the logarithm in eq. (|3.14 ) is the same as for 
the infrared-divergent piece of the corresponding one-loop resultR 


4 Gauge-dependence 


We show now that the same result (|3.10|) is obtained in general covariant gauges provided 
the large time limit in eq. (|2.25| ) is taken with an infrared cutoff in the gauge sector, in 
order to eliminate the contribution of the spurious degrees of freedom. 


In the covariant gauge of eq. (|2.11|) , the photon spectral density reads 


p(<?o,q) = 


V - <?o( v ■ q/ 

, q 2 - ql , 


Pi(qo,q ) + (i-(v-q) 2 )p t (g 0 ,g) + Ap A (go,q)- (4.1) 


The longitudinal and transverse spectral functions pi and pt are the same as in eq. (0, 
and 


Pi(«>,q) = (9o - v ' q) 2?r e(q 0 )6(q), 


(4.2) 


^Actually, a different coefficient was reported in Refs. @>0> but the difference is apparently due to 
an error in the calculations there. 
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where e(go) = 0(g o ) — 9(—q 0 ) and S'(q 2 ) is the derivative of S(q 2 ) with respect to q 2 . 

The electric and magnetic spectral functions in eq. ( |4.1| ) yield the same contributions 
to eq. ( [2. 25| ) as the corresponding functions in the Coulomb gauge (cf. eqs. ( |3T| ) and 
Q3.8| )). This is obvious for the magnetic sector. In the electric sector, the large-time limit 
introduces the delta function 5{y-q ) (see eq. ( p.271 ) )•, and the projection factor multiplying 
Pi(qo,q) ineq. (S3) is equal to one for g 0 = v ■ q. The same argument applied in the gauge 
sector seems to imply that the contribution of p\(qo,q ) does also vanish, because of the 
factor (go — v ■ q ) 2 in eq. (|4~2l). However, since the spectral function S'(q 2 ) has support 
precisely at the integration limits go = ±g, we should be more careful when taking the 
limit v ■ q —> 0. 

The contribution of the gauge sector to A (f) factorizes as exp [—A g 2 TF\(t)], where 

Fxit) = /(^/^' ,i( *>’ q) 
f d 3 g 


1 — cos t{v ■ q) 


(2tt) 


2vr g 0 
dg 0 
% 


(v ■ g) 2 

e(?o)^'(g 2 )fl-cos t(v-q) 


(4.3) 


By noting that <5'(q 2 ) = (l/2g 0 )(d5/dg 0 ), we can perform an integration by parts to 
compute the integral over g 0 . After also computing the angular integral, we obtain 




00 dg |i 


27T 2 g 2 


sin 2 qt 1 — cos 2 qt 

2 qt 4 


(4.4) 


Although this last integral is infrared finite, we nevertheless compute it with an infrared 
cutoff /i. A straightforward calculation then yields 


Fx(t) = 


8ii 2 ji 


„ „ sm 2 at 

3 — 2- — cos 2 ut 

2 [it 


t 


An 2 


si (2 /it), 


(4.5) 


where si(x) = — &z(sm.z/z) is the sine integral function [19 


If we remove the IR cutoff by letting p —> 0 at hxed t, then, by using si(0) = —vr/2, 


we get 


Fx(t) = 


8n 


so that A (t) becomes (cf. eq. (|3.10|) ) 

A(f > 1/c o p ) ~ exp |— aTt ^ln(o;pt) + 0.12652... + A/2^j. 


(4.6) 


(4.7) 


(The gauge-dependent piece in the exponent coincides with the corresponding piece of 
the one-loop damping rate, = XaT/2 [20|. 
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However, if we consider the large-time behaviour at fixed /i, then we can use the 
asymptotic expansion of si (a:), that is 

— si (a: > 1) ~ ( cosa: H-h 0(1 x 2 ) ) , (4.8) 

x \ x J 

to obtain 


F x (t > 1 /n) 


1 

87r 2 /i 



sin 2/at 
2 /it 


+ 0(l/p?t 2 ) 


(4.9) 


In this case, the sole effect of the gauge-dependent piece F\(t) at times t S> l//a is to 
change the normalization of the propagator, by a factor 

exp(— A g 2 TF x (tf) ~ exp (-A ^ ^ = z(T, fi, A), (4.10) 

which is both gauge-dependent and cutoff-dependent. 


The gauge-dependent contribution to the damping rate, eq. ( |4.6|) , arises because the 
on-shell fermion is kinematically allowed to “decay” with the emission, or the absorption, 
of a massless gauge “photon”. At T = 0, such an emission process cannot occur: by 
kinematics, the emitted photon must be colinear (6 — 0, q 0 — q), and the corresponding 
phase space vanishes. But this is not so at finite temperature, because of the Bose- 
Einstein factor N(q 0 ) ~ T/q 0 which diverges as q 0 —> 0 (see Appendix B of Ref. [Q) 
for an explicit calculation). After HTL resummation, the gauge sector is the only one 
to contain massless fields. The unphysical decay channel can be suppressed by giving 
the gauge photon a small mass /i, as originally proposed by Rebhan |2T| (see also Refs. 
). As we have seen, this procedure ensures the gauge-independence of the damping 


22 , 


rate, to the order of interest. 


Further insight may be gained by a comparison with the corresponding results at 
zero temperature 0. After ultraviolet renormalization, the retarded BN propagator at 
zero temperature is given by 


S(t, p) oc e lEpt exp j(3 — A) ln(mt)| = (mt) 1 ' 3 A) 2n e lEpt . (4-11) 

Thus, in the energy representation, the mass-shell singularity is generally a branch point, 
rather than a simple pole: 


S{p) oc 


1 

u ■ p — m 



(3-A) & 




(4.12) 


where = (uq, u) is the fermion four-velocity, tA = p M /m, with u 2 = 1. 


Note that both the physical and the gauge sectors of the photon propagator con¬ 
tribute to the mass-shell behaviour in eq. (|4.12| ): at T = 0, the gauge held quanta are 
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massless in both sectors. Furthermore, no infrared regulator is necessary: in deriving 
eqs. (|4.11|) -( |4.12| ), one encounters no IR divergences, and the position of the mass shell is 
gauge-independent, as it should. Still, if one wishes to perform soft-photon computations 
in perturbation theory in any other gauge than the Yennie gauge (A = 3), it is conve¬ 
nient to introduce an infrared regulator, so as to recover the simple-pole structure of the 
mass-shell. When the photon is given a small mass the propagator ( 4.12 ) is replaced 
by 


S(p ) oc 


u ■ p — m 


(4.13) 


where the residue 


z(fi, A) = exp {( 3 - A ) ^ ln(m/p)J , 


(4.14) 


is gauge-fixing dependent and also cutoff-dependent. It may be compared to the finite- 
temperature normalization factor in eq. ( 4. 101 ) . 


We see that, as a consequence of the Bose-Einstein enhancement of the soft photon 
processes, the divergence of the “residue” z(T,fj,, A) as p —► 0 is linear at T > 0, rather 
than just logarithmic at T = 0. Moreover, if at T — 0 the introduction of a photon mass 
is just a matter of convenience, at T > 0 the use of an infrared regulator in the gauge 
sector is compulsory in order to eliminate the contribution of the non-physical degrees of 
freedom and avoid the gauge-dependence of the mass-shell. 


5 Some results for QCD 


We consider now the generalization of the previous arguments to the non-Abelian case, 
that is, to the high-temperature, weakly-coupled (g(T) <C 1) quark-gluon plasma. The 
self-interactions of the soft gluons prevent us from getting in this case an explicit solu¬ 
tion. However, it is expected El that these interactions generate screening of the static 
magnetic fields. If this is so, the corresponding screening length, typically of order l/g 2 T 1 
provides then a natural IR cutoff of order g 2 T in the perturbation theory for 7 . 


We shall investigate this possibility in the next subsection, in the framework of a 
toy model which is QED with an infrared cutoff /j ~ g 2 T in the magnetic sector. By 
solving this model in the BN approximation, we shall obtain a qualitative picture of the 
effects of the magnetic mass on the large-time behaviour of the fermion propagator. 


Then, we shall propose a functional integral representation for the propagator of a 
hard quark or gluon which, being formulated in three-dimensional Euclidean space, is a 
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priori well suited for lattice calculations. This formulation allows for a direct numerical 
study of the particle decay in real time. 


5.1 QED with a magnetic mass 


To implement magnetic screening in QED, we replace the massless static transverse prop¬ 
agator by its massive version, 


Aj(0,q) 


&ij Qi Qj 

q 2 + /i 2 


(5.1) 


with n ~ g 2 T. Of course, such an infrared behaviour could not occur in QED, where the 
correct magnetic polarization tensor fl f (0,g) vanishes like q 2 when q —> 0, to all orders 
in perturbation theory [^]. We simply use “massive QED”, as defined by eq. (|5.1|) , as a 
crude parametrization of the non-perturbative screening effects in QCD. 


Strictly speaking, when fi > 0 we have no infrared divergences. However, as long 
as /i ^ g 2 T, the dominant contribution to the damping rate is still given by the static 
magnetic photons. For instance, to one loop order the static mode yields 

7 at = otT In — ~ g 2 T ln(l/g), (5.2) 

h 

which is enhanced by a factor ln(l/g) as compared to the contribution of the non-static 
modes (which is ~ g 2 T , as in eq. ( |3.1G| )). Moreover, the higher-loop diagrams contribute 
terms of relative order (aT//i) n_1 , where n is the number of loops, so that the perturba¬ 
tion theory breaks down for /i < g 2 T ||. This is why a non-perturbative calculation is 
necessary even in the presence of an infrared cutoff /i ^ g 2 T. 


In order to get the leading contribution to the damping factor, we can restrict 
ourselves to the interactions with static (q 0 = 0) magnetic photons J?], ||. In practice, 
such a calculation amounts to preserve only the contribution Fjji(t) in eq. (|3.6| ), where 
however the massive propagator (|5.1|) must now be used. This gives: 


A M (t) ~ exp 



d 3 g 1 1 —cosf(v-q)) 

(27t) 3 q 2 + /i 2 (v ■ q) 2 J 


(5.3) 


As explained in section 3, the integral in eq. (|5.3|) has to be computed with an upper 
cutoff u p ~ gT , to account approximately for the effect of the neglected non-static modes. 
(In eq. (|3.6Q , the upper cutoff was provided by the substracted term 1 /(g 2 + u 2 ). In 
eq. (|5.4| below, we shall End convenient to introduce this cutoff in a different way. At 
large times u p t 1, the leading contribution to the damping factor is indeed insensitive 
to the precise value of the UV cutoff, and also to the specific procedure which is used for 
its implementation H.) 
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To perform the integral in eq. (|5.3| ), we write A p (t) = exp (—g 2 TF p (t)), with 

g e *(v-q)(si-s 2 ) 


F “ w = U‘, dsi ( A ^J 


8tt Jo 

1 rt 

An J 1/, 

t 


= — / ds 


ds 2 

dsi / dsj 
Jo 

t — s 


(27t) 3 g 2 + /i 2 
g-/i|si-S 2 


Si — s 2 


0 (|si - s 2 | - 1 /^p) 


,-lis 


dx 


Q-V/Up _ Q-Ht 


An 


(5.4) 


I fi/ojp X /if J 

In this calculation, the ultraviolet cutoff has been introduced, in the second line, in the 
function 0(|si — s 2 | — l/u> p ). For the purpose of a graphical representation (see Fig. |3|), we 
rewrite the final expression above as F /t (f) = ( t/An)L(x,y ), with x = fit, y = n/u p , and 

p y _p x 


L{x,y) = Ei(y) - E^x) 


x 


(5.5) 


where E!(x) is the exponential-integral function |0|, E^x) = /“ dz (e~ xz /z). For /! ~ 
g 2 T, we have y ~ g -C 1. We recall that the above calculation only makes sense for large 
enough times, oj p t > 1 or x > ?/. 

Since the expression in eq. (|5.4|) involves two energy scales, namely u p and /i, with 
H -C oj p , we distinguish between two regimes of time: (i) very large times, i > l//i (i.e., 
x>l), where 


and (ii) intermediate times, A/u p Ct C 1/y (i.e., y - Cx<l), where 

F M - ^ (inK*) + 0(1)). 


(5.6) 


(5.7) 


Thus, the “magnetic mass” g is only felt at sufficiently large times — where A p (t) decays 
exponentially in agreement with the one-loop result ( |5.2| ) — while it has no effect at 
intermediate times. We remark at this point that, when discussing the lifetime of the 
excitation, it is rather the intermediate times which matter, since for asymptotically large 
times t > 1 /g 2 T ~ 1/g the excitation has already decayed. This behaviour, eqs. (|5.6Q — 
S), can be also observed in Fig. [|. where we have represented L(x,y ) as a function of 
x for a fixed, small, value of y (namely y = 0.01). For x of order one, one clearly sees on 
this figure the transition between the two types of behaviour, as described by eq. (|5T| ) 
and eq. O, respectively. 

It has been suggested, first by Lebedev and Srnilga [lTJ, that when computing 
the damping rate to one-loop order, the damping rate itself should be self-consistently 
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Figure 3: The function L(x, y ), eq. (|T5l) , is represented as a function of x for fixed y = 0.01 
(continuous line). We have also represented the asymptotics Li(x, y) = \n.{x/y) — l (dotted 
line) and L 2 (x, y) = ln(l jy) — ye (dashed line). These are good approximations to L(x, y) 
in the domains i/CiCl and x > 1, respectively. 
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resummed in the internal hard line. The usual argument goes as follows: since 7 ~ g 2 T is 
of the same order as the infrared cutoff /i, it should be taken into account when studying 
the infrared behaviour of the integrand. If we do that, by following Ref. (11|, then the 
one-loop result (|5.2|) is modified to 


aT 

7 ~ —— In ■ 


cut 


(5.8) 


2 / 1 2 + 2 fij ' 

(Up to appropriate color factors, the same result is obtained in QCD, for both quarks and 


gluons | 1 T[.) However, this is not correct: the self-energy resummation advocated in the 
procedure leading to eq. (|5.8|) should be accompanied by a corresponding resummation 
in the vertex function, so as to respect gauge symmetry. As discussed in Ref. ||, the 
vertex corrections generate new infrared divergences, and, when added to the self-energy 
corrections, conspire to give a leading-order estimate for the damping rate which has 
precisely the form indicated in eq. ( p. 2 |) . (See Appendix C in Ref. || for more details.) 
At this point, it might be useful to emphasize that the BN calculation provides precisely a 
self-consistent resummation of the fermion propagator near the mass-shell, together with 
the appropriate resummation of the vertex function, as required by gauge symmetry. 


5.2 QCD 

Going now to QCD, we first observe that the Bloch-Nordsieck approximation remains 
relevant to discuss the large-time (or mass-shell) behaviour of the quasiparticle propagator, 
and this for both quarks and (transverse) gluons. (We consider here a hard quasiparticle, 
with momentum p T .) While for quarks this approximation is easy to justify, by 
analogy to QED, the case of gluons requires more care and is discussed in Appendix C. 
Moreover, we expect the leading large-time behaviour to be given by the quasiparticle 
interactions with static (g 0 = 0 ) and very soft (q —> 0 ) magnetic gluons: indeed, these are 
the interactions which generate the infrared divergences of the perturbation theory ||. 

What is new with respect to QED, is that the relevant self-energy corrections also 
include the mutual interactions of the internal gluons, expected, in particular, to lead to 
magnetic screening. A typical Feynman graph contribution to the self-energy is depicted 
in Fig. [|. The continous line in this diagram is hard, and may represent either a quark, or 
a gluon. The wavy lines denote static magnetic gluons, and all the loop integrations are 
three-dimensional. As explained in section 3, these integrations involve an upper cutoff 
u; p . In QCD §, 


a 2 T 2 

c ol = {2N + N t ) 


18 


(5.9) 
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Figure 4: A generic self-energy diagram in QCD which yields infrared divergences on the 
mass-shell. The continuous line is a hard particle (quark or gluon). The wavy lines are 
soft magnetostatic gluons. All the loop integrations are three-dimensional. 


for N colours and Nf flavours of thermal quarks. (Note that the gluonic loops attached 
to soft internal lines would be finite even in the absence of the ultraviolet cutoff.) 


In the BN approximation, all such diagrams are formally resummed by the follow¬ 
ing functional integral (the gauge-fixing terms are not written here explicitly; see the 
discussion below, after eq. Q5.14Q ): 

1 


S(x — y) = Z 1 / VA G(x , y | A) exp — 


4T 




(5.10) 


where A a (x) is a static color field, F“ = diA “ — djA“ — gf abc A b Aj, ( f abc are the structure 
constants of the colour group), and G(x,y |A) satisfies the equation 


i(v ■ D lr )G(x,y |A) = a (4, (i - y), 


(5.11) 


where D p = d^ + igA^, A M = (0, A), and Aj = A“T a is a color matrix in either the adjoint 
or the fundamental representation (for gluons or quarks, respectively). The plasma effects 
do not modify the gluonic action in eq. (|5.10| ) (recall that the HTL-corrections vanish for 
static magnetic fields H), but only enter through the upper cutoff u p gT. 


The solution of the BN equation ( |5.11|) with retarded boundary conditions is im¬ 
mediate [0: 


Gr{x,v | A) = -i6(x 0 - y 0 )S {3) (x - y - v(x 0 - y 0 ))U(x,y) 
f/(x, x — vt) = P exp jig J ds v ■ A(x — v(£ — s))j , (5.12) 

where the path-ordering operator P is necessary since the color matrices A(x) at different 
points along the path do not commute with each other. 


The retarded propagator Sr(x — y) is calculated by inserting eq. (|5.12|) in the func¬ 
tional integral ( |5.10| ). It can be written as 


S a R b (t, p) = —i 5 ab 9(t) e _i ^ v ' A(£), 


(5.13) 
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with 


A (() = A f J PA Tr f/(x, x — vi) exp j — /— J d 3 x B 


(5.14) 


B 2 = BfBf, and = (1/2 )e i jkF°' k is the chromomagnetic field. In order for eq. (|5.14|) to 
be well-defined, it is further necessary to choose a gauge within the functional integral. 
(Recall that the parallel transporter Tr f/(x, x — vt) is not invariant under the gauge 
transformations of the background color field.) However, we shall argue below that the 
dominant large-time behaviour should be independent of the gauge-fixing condition. 

The functional integral (|5.14|) cannot be computed analytically because B a is non¬ 
linear in the gauge potentials. However, we may expect the large-time behaviour of A (t) 
to be similar to that of the model discussed in section 5.1, that is (cf. eqs. (|5.6|) -- (|5.7|) ), 

g 2 T 


A(l/u ) p « t <C 1/fj,) ~ exp <j — C r —— t (ln(o;pi) + 0(1)) 


(5.15) 


at intermediate times, and, respectively, 

A(i > l//i) ~ exp < — C r 


g^T 

An 


t (in + 0(1)) 


(5.16) 


at very large times. In these equations, C r is the Casimir factor of the appropriate 
color representation (i.e., C q = (N 2 — l)/2 N for a hard quark, and C g — N for a hard 
gluon), and the magnetic mass /i ~ g 2 T is expected to come out from the soft gluon 
mutual interactions. Note that the 0(1) terms in the above equations are not consistently 
determined by the present approximation: Indeed, from the experience with QED, and 
also from the one-loop calculations in QCD E3 |TT| , p~2| , we know that such terms receive 
contributions from the non-static gluon modes, and that they may be sensitive to the 
gauge-fixing condition (cf. section 4). 

To verify this picture, one could rely on a lattice computation of the Euclidean 
functional integral (|5.14|) . The parallel transporter Tr f/(x, x — vt) is easily implemented 
as a product of link operators. (Recall that v is a fixed unit vector, e.g., v = (0, 0, 1), so 
that H(x, x — vt) is a product of link operators in the z direction, from z 0 — t to z 0 , with 
Zo an arbitrary site on the lattice.) Since the expression (|5.14j) is defined with an upper 
cutoff c o p ~ gT , the lattice spacing a is fixed: a ~ l/cUp. From perturbation theory, we 
expect the decay of A (t) at times t a to be only logarithmically sensitive to the precise 
value of a (cf. eqs. (|5.15|) -( |5.16|) ). 

The objective of a lattice calculation would be then to verify the large-time be¬ 
haviour predicted in eqs. ( |5.15| )- (|5.16|) ). By observing the interplay between these two 
regimes, one may verify what is the typical scale for the emergence of magnetic screenig: 
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indeed, we expect the transition between two regimes to occur for t ~ 1/p. However, this 
conld not be sufficient for a quantitative measure of the magnetic mass, because of the 
theoretical uncertainty on the subleading term 0(1). 

The main limitation against an explicit calculation comes from the lattice size: 
indeed, in order to verify the aforementioned picture, one needs a small coupling constant 
g <C 1 — to ensure a clean separation between the scales g 2 T and gT — together with 
large values of time, up to t max 1 /g 2 T ~ ajg. Thus, the lattice should have at least 
N 3 sites, with N t max /a 1/g. 


6 Conclusions 


In this paper, we have completed the analysis of the Bloch-Nordsieck propagator in hot 
QED, and we have also discussed the usefulness of such an approximation for a high- 
temperature QCD plasma. 


As compared to Refs. am several points have been clarified by the present anal¬ 
ysis. First, the three-dimensional model of [7|. | 8 [] suffers from a spurious ultraviolet diver¬ 
gence coming from the restriction to the static photon mode. We have shown here that 
the contribution of the non-static modes provides a dynamical cutoff at momenta rs./ gT. 
Not only this justifies the cut-off procedure used in Refs. 11 , but it also allows one to 
compute explicitly the subleading term in the large time behaviour (cf. eq. ( |3.10| )). 

Secondly, the effects of the gauge-fixing procedure only enter at the level of the 
subleading term. Thus, by computing this term in different gauges, one can study the 
gauge-(in)dependence of the large-time decay. We have performed this computation in the 
Coulomb gauge, and in a generic covariant gauge, with conclusions which agree with the 
one-loop calculations in Refs. [^0[ [2lJ; namely, the subleading term is gauge-independent 
if computed in the presence of an infrared cutoff in the gauge sector. Physically, such a 
cutoff separates the particle mass-shell from the threshold for the spurious emission or 
absorbtion of massless gauge photons. With a non-zero infrared regulator /i, the gauge- 
dependent contribution to the damping vanishes in a vicinity ~ p of the mass-shell, or, 
equivalently, for times t 1/p. Thus, to avoid spurious gauge contributions over a 
particle lifetime, one should choose p ^ 7 , where 7 p 2 Tln(l/ g). 


Concerning QCD, we have argued that the BN model may still be a relevant ap¬ 
proximation for the study of the quasiparticle mass-shell. As compared to QED, the 
non-Abclian model is complicated by the soft gluon self-interactions. The problem sim¬ 
plifies considerably when one considers only the dominant contribution due to the static 
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magnetic gluon modes: then, not only all the HTL corrections vanish, but the needed 
path-integral can in principle be computed on a three-dimensional lattice, with a hxed 
lattice spacing a rs./ i /gT. 
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Appendix A 


In this appendix, we construct the BN propagator for a thermal fermion with mo¬ 
mentum p ~ T. The final result turns out to be essentially the same as that obtained in 
section 3 in the case p> T. 


We use the imaginary-time formalism which has been developed in Ref. || . In this 
formalism, we have to solve the imaginary-time BN equation (cf. eq. ( |2.10|) ): 


- 0 • D x )G(x,y\A) = S E (x, y ), 


(A.l) 


with antiperiodic boundary conditions (cf. eq. (|2.14|) ): 

G e (t x = 0, T y \A) = - G e (t x = (3, T y \A), (A.2) 


and similarly for r y . In these equations, the time variables are purely imaginary (xq = 
t 0 - ir x and y 0 = t 0 - ir y , with 0 < t x , r y < (3 and S E (x - y) = S(t x - r y )h(x - y)), and 
the gauge fields are periodic in imaginary time: A p {r = 0) = A p {r = (3). In Ref. ||, 
we have solved this equation explicitly in perturbation theory, i.e., as a series in powers 
of gA, and then we have performed the functional integration over the gauge holds (cf. 
eq. (|2.9|) ). The resulting propagator can be written as (cf. eq. (|2T5|) ) 


iS(x 0 -y 0 ,p) = 9 (t)S > (t,p) - 9(-t)S < (t,p), (A.3) 

where the analytic functions S > and S < are obtained in the form Q 

S' < (r, p) = e~ rEp V ( E p ; u = — r) for — [3 < r < 0, 

S > (t, p) = e {0 ~ T)Ep V(E p ]u = (3 - t) for 0 < r < (3. (A.4) 

In this equation, E p = v ■ p is the BN mass-shell and the function V(E p ; u) is given, for 
0 < u < f3, as a formal series in powers of g 2 [|| : 


V(E p ,u) 


'«( v ■ P) + /[dgidg 2 ••• d q n ] 

n> 1 ^‘ 


^(gl)-P(g2) --P(gn) 

(v ■ q\) 2 (v ■ q 2 ) 2 ... (v ■ q n ) 2 
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n(v ■ p) — n(v ■ (p + qi)) e u ( v ' qi ' ) — n(v ■ (p + q 2 )) e u ( v ' q2 '> + 
... + (-l) n n(v ■ (p + qi + q 2 + - + q n )) e -«"-(?i+®>+~+*Ol , 


(A.5) 


with D(q) = v M *D tlu (iu Tn , q)v u . In this and the following equations, the photon energies 
g 4 ° are discrete and purely imaginary: q° = iuj m = i2nmT, with integer m (Matsubara 
frequencies). The measure in the momentum integrals is denoted by 

(A.6) 


The thermal factors make the momentum integrals in eq. ( A.5 ), like 

/(Cy n(v ' (p+q))e " <vq) ' 


(A.7) 


convergent for any 0 < u < (3. This, in turn, ensures the analyticity of the functions S < (r) 
and S > (t) in ( |A.4j) jl6[| . By analytically continuing these functions toward the real-time 
axis (i.e., by replacing t —* it with real t), one constructs the retarded propagator: 

Sn(t, p) = —i6(t)(S > (t, p) + S < (t, p)j 

= -iO(t)e- itEp {e 0Ep V(E p ;u = [3-tt)+ V(E p -,u =-it)}. (A.8) 

Note, however, that the analytic continuation to real time can be done only after per¬ 
forming the Matsubara sums in all the terms of the series in eq. O- 

Fortunately, this can be done easily in the relevant regime of large-time (t 1/gT). 
According to the discussion in section 3, we expect then the momentum integrals to be 
dominated by .soft photon momenta, q ^ gT. Indeed, the photon propagator D(q), which 
can be rewritten as (with p(q) = v fl p^(q)v l/ ; cf. eq. (|Q|) ) : 

r°° do; p(u, q) 


D(iu m , q) = 


(A.9) 


l—oo 27 t uj — ioj m 

provides, through the spectral density p(u;,q), an effective upper cutoff ~ gT for the 
integrals over q. (Recall that the functions fli(u,q ) and f3 t (uj,q) in eq. ( fl.2|) are rapidly 
decreasing for q gT.) Strictly speaking, this cutoff becomes effective only after u is 
continued to (3 — it or —it. However, we may anticipate for its effect and supply the 
integrals over q^ in eq. ( |A.5| ) with an upper cutoff ~ gT. Then the photon momenta are 
limited to values |q| |p| ~ T, and we can replace n(v • (p + q)) by n(v • p) up to terms 

of order q/T ^ g. The fermion occupation factor n(v ■ p) in eq. (|A.5| ) then factorizes, 
and the resulting expression can be resummed into an exponential: 


V{E p -u) 


-‘vi 

A<| 


n(E p ) A < (u) 


u 


= exp 


D(q) 


X _ e -“h-9) 
(y ■ q) 2 


(A.10) 
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At this stage, we can then perform the Matsubara sum over q° = ico m (by using the 
spectral representation in eq. (|A.9|), together with contour integration), and obtain 


A < (u) = exp < — g 2 


d 4 q 

( 2 W 


p{q o, q) 


1 _ p -u{y-q) v 

(1 + N(q 0 )) * ^ - (1 + N(y ■ q)) 


{v ■ qf 


v ■ q 


(A.ll) 


where v ■ q = go — v ■ q (we have renamed go the real energy u). The last expression can 
be now continued to u —» — it , with the result: 


A < (—it) = exp < — g 


d 4 q 


K - (2tt) 4 
= exp{— A(t), 


p{q) 


1 _ P it{v-q) Af 

(1 + N(qo)) — -— + (1 + N(y ■ q)) 


(v ■ qy 


v ■ q 
(A.12) 


which involves the same functions $(t) and A(t) as in section 3 (cf. eq. (|2.24|) and ( p.25| )). 
At this point, the momentum integral in eq. (|A.12|) is ultraviolet hnite and the cutoff can 
be removed. Recalling eq. (|A-4j) , we can hnally write: 


5 ,< (t,p) = e ltEp n(E p ) A(t), 

where we have ignored the inconsistent phase $(i). 

To compute S > (t, p), we start with (cf. the second eq. (|A.4|) ): 


(A.13) 




V(Ep, u = /3 - t) « [1 - n(E p )\ A > (r) 


A>(r) = exp -g- 


D(q) 


*r(v-q) 


(v ■ qf 


(A. 14) 


where 1 — n(E p ) = e l3Ep n(E p ) has been factorized by the same approximations as above. 
After performing the Matsubara sum and the analytic continuation r —> it , we finally 
obtain (within the same accuracy as in eq. ( |A.13|) ): 


S > (t, p) = e~ ltEp [l — n(E p )\ A(t). 


(A.15) 


Thus, for sufficiently large times, both functions S < (t ) and S > (t) decay like A (t), eq. (|2.25| ). 
The same is therefore true for the retarded propagator, as given by eqs. ( |A.8| ), (|A.13| ) and 

(CT): 


l<Sfl(t,p)| OC A(t), (A.16) 

which is the result quoted in eq. ( p.33|) . 
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Appendix B 


In this appendix, we calculate the double integral in eq. ( |3.7| ) , thus proving the 
result quoted in eq. (|3.8|). The method to be used here was suggested to us by Jean-Yves 
Ollitrault (see also EBU). We first write 


/ = 


dqq 


rq dq 0 


% 2 


-q 27Tgo 


A(<7 o,Q) — -y AAoW) + Vt(qo, O') = h + h + hi (B.l) 


T 


where the three pieces I s , s = 1, 2, 3, correspond to the three terms within the integrand. 
To illustrate the method, we compute the second piece in detail: 

r q dq 0 ql 


h = - 


dqq 

Ay 

q l 


-q 2vrg 0 q‘ 


■Pt(qo,q) 


^QoQo^Qo -Q 2 - n t (go,?)j 


In going to the second line, we have used the familiar sum-rule |16 

/ % PtiMih q) = 1 , 

J — oo Zll 

together with the parity property Pt(—qo , o) = ~Pt(q o, q), to write 


I- A (3o, 3) = 1 -2jf ^9oPi(«o.9); 


d«o 


(B.2) 


(B.3) 


(B.4) 


then we have related the on-shell magnetic spectral density to the plasmon pole in the 
transverse photon propagator: p t (qo > q) = 27rJ(gQ — q 2 — TI t {q 0 ,q)). We also recall that, 
in the hard thermal loop approximation, n t (g 0 , q) is a function of q 0 /q alone. 


The integral over q in eq. ( |B.2| ) is well dehned as it stands. However, in order 
to work out separately the two terms within the braces, it is necessary to introduce, at 
intermediate steps, an ultraviolet cutoff A and also an infrared cutoff //. The first term 
reads then 


MA.p) = - / A - = - In-. 

Jn q p 


The second term, 

/ 22 (A,aO = 2 f — f dqoq 0 $(qo - q 2 -~R t {qo 

Ju a Ja v 


(B.5) 


(B-6) 


involves an integral along the transverse plasmon dispersion relation, q 0 = uj t (q) with 
to 2 (q) = q 2 + n t (ujt/q). To perform the integral, we use the following change of variables: 

da: dy 


x = q 0 /q, 


y = q o - 9 2 , 


dqdq 0 = 


2(x 2 - 1 ) ’ 


(B.7) 









and get 

/ , r x r x M 

dx dy— — -6(y-n t (x))= dx 

J X i. J Xm. 


X 


= - In 


%M ~ 1 


X z 


r- — 1 
1 


(B.8) 


The integration limits x m and Xm are obtained as follows: As q —> /x, x —> u t (n)/H- The 
dispersion relation ag(g) can be found, e.g., in Refs. [18|. 1| [RJ. For y —> 0, uy^/i) —> u p , 
and x —> xm(ix) = u; p //n As q —> A (with large A u; p ), u; 2 (A) ~ A 2 + 3a; 2 /2, and 


£ —> x m (A) = 1 + 3tn 2 /4A 2 . Together with eqs. (|B.5|) and 


, this gives 


r r r , A 1 2A 2 12 

J 2 = / 2 i + / 22 = - In - + - In — = - In - . 

H 2 3/R 2 3 


(B.9) 


The remaining integrals I\ and I 3 are evaluated similarly. In the process, we need 
the following sum rules B0 


dq ° Pi(qo,q) = 4 - — , 4 . .9 - 2 / + 


l-q 2vrg 0 
f q dq 0 
l-q 2tt g 0 


g 2 g 2 + 3a; 2 


<?o 




/°° dgp / 2 2 t-t / 

/ - o[q 0 -q - n t g t 

Jq qo v 


<? 2 + ^ 2 


(B.10) 


where 11/ and 11, are the polarisation functions in the hard thermal loop approximation. 
(We use the same notations as in Ref. ||.) In the computation of I 3 — which involves 
u t (q 0 , q ) — we change the integration variables as in eq. (El above, and obtain: 


h = - In - 


2 ’ 

which happens to cancel J 2 , eq. As for the electric piece I \, we write 


(B.ll) 


h = 


dgg< -y 


2/ ^(g 2 + n z (g c 


g 


= In 






g 2 + 3a; 2 
dg 0 


qo 


A 


— 2 f°° dqq f°° S(q 2 + Hi(qt 
•'m J q g 0 v 


(B.12) 


where an infrared cutoff /1 was introduced when separating the terms inside the braces. 
In the second term, we change the variables according to 


and get 


x = qo/q, 


- f — [ dyS(y + Ui(x)) = -lnx M = - In ^ . 
Ji x J n 


y = q 2 , 


dqdq 0 = -dxd y, 


(B.13) 


(B.14) 


The upper limit was obtained as %(/i) = og(/i)//i ~ u p /y for fi —> 0. From eqs. (E3) 
and (|B.14|) , we finally obtain 


h = 2 ln3 > 


(B.15) 
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which is the result quoted in eq. (0). 

We finally evaluate the momentum integral in eq. ( |3.12| ) in the non-relativistic limit 
u€l. Since |go| < vq <C g, we need the spectral functions /?z,<?) only for very small 
frequencies JTB| : 


A (go < g) - 
A (go < g) - 


3t tuj 2 p (g 0 /g) 

(g 2 + 3 lu 2 ) 2 
3?r^p (go/2g) 
g 4 + (37ru;2g 0 /4g) 2 ' 


(B.16) 


Corresponding to the three terms in eq. ( |3.12| ), we write C{v) = C\(y) + C 2 (v) + C 3 (v). 
The electric contribution is evaluated as follows: 


Ci(v) = 


3 , [ vq dg 0 
2 i d gg / o— 

v z Jo J-vq 2nqo 


at ^ 1 /'°°u 

A(go,g) --n dgg 

U 2 JO 


(g 2 + 3 oj 2 ) 2 2u 


= 7T> (B- 17 ) 


where the neglected terms are smaller, at least, by two powers of v (since A(go,g) is an 
odd function of g 0 ). 

The first magnetic contribution is 


Co(v) = 


1 r°° ri dg 0 ql R , s 
— dgg / -- 2 A(go,g) 

u 2 Jo J-vq 2 tt g 0 g 2 


?>Tt 2 U 2 V JO 


dgg 


1 - 


arctan y(g; r;) 
y(g;u) 


(B.18) 


where we have used the approximate expression ( |B.16|) for A (go, g) to perform the integral 
over g 0 , and we have denoted ?/(g; u) = 37rcj 2 u/(4g 2 ). In the remaining integral over g, we 
make the obvious change of variables y(q',v) = t, with dq/q = dt/(2t), and obtain 


Co(v) = 


r°° d t 
t 2 


arctan t 


t 


1 

7T JO 


7r Jo 

which is independent of u. 

Finally, the second magnetic contribution reads 


roc 1 

J o dy [1 - yaiccoty] = -, (B.19) 


CsW = / dgg f ^-^(g 0 ,g) - [ 

JO J-vo 2-KQn JO 


-«« 27r go 


dg 

g 


7T 


arctan y(q',v) — 


cut 


g 2 +(n 2 


(B.20) 


where we have also used the definition flOl) of u t (q 0 ,q). Note that, for any v > 0, the 
remaining integral over q is well-defined, and saturated by soft momenta, q ^ cu p . Still, 
the limit v —> 0 is not well-defined (because of potential infrared singularities), so that we 
need to perform the momentum integral before studying the small v behaviour. By using 
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(a) 




Figure 5: Some two-loop self-energy corrections in QCD. The continuous line represents 
a hard gluon, which is nearly on-shell. The wavy lines are soft magnetostatic gluons. In 
the on-shell limit, diagram (b), which also involves a four-gluon vertex, is less infrared 
singular then diagram (a). 


the same change of variables as above, we rewrite eq. (|B.2C| ) as 

C 3 (v) ~ f 
Jo 


r°° dt 

\2 

t ] 

Jo 21 

— arctant — 
.7r 

t + V. 


, arc tan t 

/ d t - 

/ 0 TTt 


dt 

IT' 


arctant — — 

h lit 

L 2j 


+ - In®, 


(B.21) 


with v = 37tu/ 4. The two integrals in the second line mutually cancel, as can be seen by 
changing t —> 1/t in any of them, and then using arctanl /t = 7t/ 2 — arctanf. Finally, 


C 3 (v) ~ 


- In v. 
2 


(B.22) 


By putting together the above results in eqs. 
result quoted in eq. (|3.13|) . 


>.171 ), ( |B.19| ) and (|B.22|) , one obtains the 


Appendix C 

In section 5 above, we have used a non-Abelian version of the Bloch-Nordsieck model 
to study the interactions between hard quasiparticles (quarks or gluons) and soft virtual 
gluons in hot QCD. In this appendix, we examine the validity of this approximation for 
the case where the hard quasiparticle is a transverse gluon. 

The dominant contributions to the hard (p ^ T) gluon propagator near the mass- 
shell at u = p (i.e., the leading infrared divergences for u —> p) come from the diagrams 
illustrated in Fig. £|. The continous line there represents the hard gluon and the wavy lines 
denote very soft (q <C gT ) static ( q 0 = 0) magnetic gluons. Once again, our strategy is to 
consider first the interactions with a classical, static, color field A a (q). Then, a typical 
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co, p co, p+q x 


Figure 6 : A typical diagram contributing to G(x, y | A) to order g n in perturbation theory. 
This diagram involve n gluon field insertions, and n + 1 free propagators Go (including 
the external lines). The external fields are purely static and magnetic. 


diagram looks like in Fig. | 6 |. The self-energy corrections in Fig. (| will be eventually 
recovered by functional integration, as shown in eq. fl5.10| ). 

A noteworthy feature of Figs. | 6 | and Q is that the hard particle is involved only in 
three-gluon (but not in four-gluon) vertices. (Of course, the four-gluon vertices do also 
enter the self-energy diagrams — see, e.g., Fig. [4] — , but they couple only soft internal 
gluons; cf. eq. ( |5.10| ) .) The reason is that, to a given order in perturbation theory, the 
diagrams which involve the hard particle in four-gluon vertices are less infrared singular. 
This can be easily verified by power counting: Consider, e.g., the two two-loop graphs 
in Fig. ||. For uj — p, the diagram in Fig. ^|a, with only three-gluon vertices, gives rise 
to a linear infrared singularity. That is, its contribution to the damping rate is of the 
order 7 ^ ~ g A T 2 /(i (up to logarithms of gT // 1 ), which for p ~ g 2 T gives 7 ^ 2a > ~ g 2 T ; 
i.e., it is of the same order as the one-loop contribution. (This leading divergence can be 
isolated by using the simplified BN Feynman rules to be derived below. See Appendix 
C in Ref. || for a detailed analysis.) The diagram in Fig. |5].b, which also involves one 
four-gluon vertex, may give rise, at most, to logarithmic mass-shell singularities. We thus 
expect 7 ^ ~ g A T 2 /p ~ g A T , which stands beyond our present accuracy, and should be 
discarded for consistency. We shall verify shortly that, for the problem at hand, neglecting 
the four-gluon vertices is indeed consistent with gauge symmetry. 

Consider the diagram in Fig. with only three-gluon vertices. The latter are linear 
in the external gluon momenta: 


- igf abc Tiji{ P, q, k) = -igf abc ((p - q)^ + (q- k)i5ji + (k - p)j8 u ), (C.l) 

where p + q + k = 0. Remember that all the external lines in Fig. [] are of the magnetic 
type, so that we need just the spacial components of the vertex function. Furthermore, 
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color indices play no special role for the subsequent kinematic approximations, and will 
be omitted in intermediate formulae. 


For all the vertices in Fig. one of the external momenta is soft, since it is carried 
by the classical color field. If q is the soft momentum in eq. ( |C.1|) , then 

Fiji(p,q,k) - r ij7 (p,0,-p) = piSij + PiSji - 2p j s u . (C.2) 

Since the approximate three-gluon vertex (|C. 2 | ) is independent of the soft momentum q, 
the Ward identities are consistent with setting the four-gluon vertex to zero, which is 
what we did before. 


Consider now a typical internal gluon line in Fig. j 6 j: It is necessarily hard and nearly 
on-shell. In the Coulomb gauge, the associated propagator reads (recall that qo = 0) 


D i:j {uj, p + q) 


&ij (pi “1“ Qi)(.Pj T Qj) ^ ^ij ViVj 

to 2 - (p + q) 2 2p UJ - v ■ (p + q) ’ 


(C.3) 


where pi = pi/p = u* and the approximate equality holds since g€j) and ~ p. That is, 
Dij(u>, p + q) ~ ( 1/2 p)Vij Gq(uj, p + q), where P l3 = S i3 — v^j is a transverse projector, 
and Go is the BN propagator (cf. eq. 0)) : 


G 0 (o;,p + q) 


1 

uj — v • (p + q) ' 


(C.4) 


In Fig. | 6 ;. all the three-gluon vertices like ( |C.2| ) appear between projectors like V l3 . By 
using the identity 


'Pim r'm i jre(P, 0, p) Pnl 2 (Sn VjV[ )Pj 2p Pn Vj , (C.5) 

it can then be easily verified that the leading contribution of the diagram ^ to the hard 
gluon propagator can be evaluated with the following simplified Feynman rules (we rein¬ 
troduce here the color indices): (i) the hard particle propagator 5 ab G 0 (u, p + q), and (ii) 
the hard particle-soft gluon vertex igf abc Vi. These are the Feynman rules which have been 
used to define the non-Abelian Bloch-Nordsieck model in section 6 . For a hard quark, the 
color indices in the above Feyman rules should be replaced by the corresponding indices 
in the fundamental representation. 
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